Slow dynamics in gelation phenomena: Prom chemical gels to colloidal glasses 
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We here discuss the results of 3d MonteCarlo simulations of a minimal lattice model for gelling 
systems. We focus on the dynamics investigated by means of the time autocorrelation function 
of the density fluctuations and the particle mean square displacement. We start from the case of 
chemical gelation, i.e. with permanent bonds, and characterize the critical dynamics as determined 
by the formation of the percolating cluster, as actually observed in polymer gels. By opportunely 
introducing a finite bond lifetime Tb, the dynamics displays relevant changes and eventually the 
onset of a glassy regime. This has been interpreted in terms of a crossover to dynamics more typical 
of colloidal systems and a novel connection between classical gelation and recent results on colloidal 
systems is suggested. By systematically comparing the results in the case of permanent bonds to 
finite bond lifetime, the crossover and the glassy regime can be understood in terms of effective 
clusters. 
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I. INTRODUCTION 

The gelation transition transforms a viscous liquid into 
an elastic disordered solid. In general this is due to the 
formation in the liquid phase of a spanning structure, 
which makes the system able to bear stresses. 
In polymer systems, this is due to chemical bonding, that 
can be induced in different ways P, , producing a poly- 
merizations process. As firstly recognized by Flory, the 
change in the viscoelastic properties is directly related 
to the constitution inside the sol of a macroscopic poly- 
meric structure, that characterizes the gel phase. In ex- 
periments Q the viscosity coefficient grows as a power 
law as function of the relative difference from the critical 
polymer concentration with a critical exponent k. The 
onset of the elastic response in the system, as function 
of the same control parameter, displays a power law in- 
creasing of the elastic modulus with a critical exponent /. 
As implicitly suggested in the work of Flory and Stock- 
mayer |l| , the percolation model is considered as the basic 
model for the chemical gelation transition and the macro- 
molecular stress-bearing structure in these systems is a 
percolating network 0, Q, • In experiments the gelling 
solution typically displays slow dynamics: the relaxation 
functions present a long time stretched exponential de- 
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cay ~ e ^ T » ' as the gelation threshold is approached. 
In particular at the gel point the relaxation process be- 
comes critically slow, and the onset of a power law decay 
is observed p. 

In many other physical systems where aggregation pro- 
cesses and structure formation take place, gelation phe- 
nomena can be observed. Typically, these are colloidal 
systems, i.e. suspensions of mesoscopic particles inter- 
acting via short range attraction. These systems are 
intensively investigated due to their relevance in many 
research fields (from proteins studies to food industry). 



Due to the possibility in experiments of opportunely tun- 
ing the features of the interactions, they also play the role 
of model systems. 

In these systems strong attraction gives rise to a diffu- 
sion limited cluster-cluster aggregation process and may 
produce gel formation (colloidal gelation) at very low 
density as a spanning structure is formed Q- The latter 
is generally quite different from the polymer gels case @ 
however, for what the viscoelastic behavior is concerned, 
this gelation transition closely resembles the chemical 
one, observed in polymer systems 9]. With a weaker at- 
traction at higher densities a gelation characterized by a 
glass-like kinetic arrest |0, ma y t> e observed. The 
relaxation patterns closely recall the ones observed in 
glassy systems and are well fitted by the mode-coupling 
theory predictions for super-cooled liquids approaching 
the glass transition On the theoretical side the ap- 
plication of the mode-coupling theory to s yste ms with 
short range attractive interactions |13L [lj, (attrac- 
tive glasses) has been recently considered and the con- 
nection with the colloidal glass transition has been pro- 
posed. The short range attraction enhances the caging 
mechanism characteristic of glassy regimes in hard sphere 
systems and produces a glassy behavior at lower densi- 
ties, depending on the attraction strength. 

At lower densities, the role of the structure formation, 
as directly observed in some systems |l6j . might be rel- 
evant in the dynamics but it has not been clarified yet. 
Also the eventual crossover to the glassy regime 0, as 
the density is increased, is not completely understood. In 
this paper we investigate the connection among colloidal 
gelation, colloidal glass transition and chemical gelation. 
Some preliminary studies have been reported in |17|. 

We have considered a minimal model for gelling sys- 
tems and performed extensive numerical simulations on 
3d cubic lattices. In Sect. [H]we give the details of the 
model and the simulations. In Sect.s ITU and ITVl the re- 
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suits relative to relaxation and diffusion properties are 
presented and discussed. In the last section some con- 
cluding remarks are given. 



II. DESCRIPTION OF THE MODEL 

A. Permanent bonds 

Our model consists in a solution of monomers diffus- 
ing on a cubic lattice. As in most experimental polymer 
gels, we choose the monomers to be tetrafunctional. Each 
monomer occupies a lattice elementary cell, and therefore 
eight vertices on the lattice. To take into account the 
excluded volume interaction, two monomers cannot oc- 
cupy nearest neighbors and next nearest neighbors cells 
on the lattice, i.e. nearest neighbor monomers cannot 
have common sites. At t = we fix the fraction cf> of 
present monomers with respect to the maximum num- 
ber allowed on the lattice, and randomly quench bonds 
between them. This actually corresponds to the typical 
chemical gelation process that can be obtained by irra- 
diating the monomeric solution. We form at most four 
bonds per monomer, randomly selected along lattice di- 
rections and between monomers that are nearest neigh- 
bors and next nearest neighbors (namely bond lengths 
1 = 2, 3). Once formed, the bonds are permanent. 

For each value of (j> there is an average value Nt,(4>) 
of the fraction of formed bonds with respect to all the 
possible ones, obtained by averaging over different initial 
configurations. 

Varying (j) the system presents a percolation transition 
at (f) c = 0.718 ± 0.005 [Tg. The critical exponents found 
at the transition agree with the random percolation pre- 
dictions (e.g. for the mean cluster size 7 ~ 1.8 ±0.05 
and for the correlation length v ~ 0.89 ± 0.01 in 3d [l8|). 

The monomers diffuse on the lattice via random lo- 
cal movements and the bond length may vary but not 
be larger than Iq according to bond-fluctuation dynam- 
ics (BFD) [23, where the value of Iq is determined by 
the self-avoiding walk condition. On the cubic lattice we 
have Iq = \/T0 in lattice spacing units and the allowed 
bond lengths are I = 2, VE, 3, VlO Hlj. We let the 
monomers diffuse to reach the stationary state and then 
study the system for different values of the monomer con- 
centration. 

This lattice model with permanent bonds has been in- 
troduced to study the critical behavior of the viscoelastic 
properties at the gelation transition|22|. The relaxation 
time is found to diverge at the percolation threshold cf> c 
with a power law behavior [Tsl ]. The elastic response in 
the gel phase has been studied by means of the fluctu- 
ations in the free energy a nd g oes to zero at 4> c with a 
power law behavior as well ,23j. 



B. Bonds with finite lifetime 

Colloidal gelation is due to a short range attraction and 
in general the particles are not permanently bonded. To 
take into account this crucial feature we introduce a novel 
ingredient in the previous model by considering a finite 
bond lifetime Tf, and study the effect on the dynamics. 

The features of this model with finite Tf, can be real- 
ized in a microscopic model: a solution of monomers in- 
teracting via an attraction of strength —E and excluded 
volume repulsion. Due to monomer diffusion the aggrega- 
tion process eventually takes place. The finite bond life- 
time Tf, is related to the attractive interaction of strength 
-E, as n ~ e E ' KT . 

We start with the same configurations of the previous 
case, with a fixed <f> where the bonds have been randomly 
quenched as described above. During the monomer diffu- 
sion with BFD at every time step we attempt to break 
each bond with a frequency 1/rf,. Between monomers 
separated by a distance less than Iq bonds are then 
formed with a frequency /f, j24[. In order to obtain 
monomers configurations that are similar to the ones with 
permanent bonds, for each value of T{, we fix /f, so that 
the fraction of present bonds coincides with its average 
value in the case of permanent bonds, Nb((j)) |2(j. 

With respect to the case of permanent bonds we no- 
tice that, as the finite bond lifetime Tf, corresponds to an 
attractive interaction of range Iq, it actually introduces 
a correlation in the bond formation and may eventually 
lead to a phase separation between a low density and a 
high density phase: There is no evidence of phase separa- 
tion for the values of Tf, and /f, considered in this paper. 
This is evident in Fig^ where typical equilibrium con- 
figurations with <fi = 0.6 are shown in two different cases: 
in Fig^a) we have the case considered in this paper, 
obtained with Tf, = 100 and /f, = 0.02; in Figm&) with 
Tf, = 2 and /f, = 1 the phase separation may occur. The 
choice of monomers of functionality 4, also in this case 
of finite bond lifetime, may correspond to a directional 
effect in the interaction |27|. 

The case of extremely large Tf, gives rise to different sit- 
uations depending on the initial condition and the bond 
creation process. We consider two extreme cases: I) 
Start with the initial configuration where the monomers 
are randomly distributed, and the bonds are randomly 
quenched. At later times the frequency of forming bonds 
is zero. This case corresponds to the permanent bond 
case (chemical gelation), described in section 111 Al II) 
Start with a random configuration of monomers and let 
the monomers diffuse and form bonds with infinite life- 
time and frequency /f,. This phenomenon of irreversible 
aggregation, with the occurrence of gelation after a span- 
ning cluster is formed, corresponds to cluster-cluster ag- 
gregation class of models for ft, = 1 [28| , with a tendency 

towards cluster-cluster reaction limited aggregation pro- 

1 1 ... 
cess [29] for < 1. This out of equilibrium phenomenon 

can be representative of colloidal gelation and will not 

be considered here. In chemical gelation and colloidal 
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gelation the formation of the critical cluster should pro- 
duce the slow dynamics. The main difference is due to 
the fact that the critical density and the temperature are 
much lower in colloidal gelation than in chemical gela- 
tion, moreover the fractal dimension is related to cluster- 
cluster aggregation models and not to random percola- 
tion. 



III. RELAXATION FUNCTIONS 

In order to investigate the dynamic behavior we study, 
for both the permanent bond and the finite bond lifetime 
cases, the equilibrium density fluctuation autocorrelation 
functions, f$(t), given by 



IS) 



< Pq -(t + t')p_ ? (t') > 

< \pS')\ > 2 



(i) 



where p<f(t) = J2iLi e~ l ^' Ti( - t \ ^(t) is the position of the 
i — th monomer at time t, N is the number of monomers 
and the average (...) is performed over the time t' . Due 
to periodic boundary conditions the values of the wave 
vector q on the cubic lattice are q — ^■(n x ,n y7 n z ) with 
n x ,n y ,n z — 1, L/2 integer values (in our simulations 
a cubic lattice of size L = 16 has been considered) [3(j ■ In 
the following we discuss the data obtained for q ~ 1.36 
(q = (it/4, 7r/4, 7t/4 )). Qualitatively analogous behav- 
iors have been observed for (q — (ir/2, n/2, tt/2 ) and 
(-7T, 7r,7r)). Undoubtly, due to structure formation over 
different length scales, a detailed study of the geometric 
properties and the dynamics for different wave vectors 
might be relevant [3l|. 

In the case of permanent bonds the system is consid- 
ered at equilibrium when both the diffusion coefficient 
has reached its asymptotic limit, and the autocorrelation 
functions have gone to zero. For <fi < <fi c we are always 
able to thermalize the system, instead for cj> > <p c it re- 
mains out of equilibrium, and it is possible that it is in 
an aging regime [3l| . In FigEl fq{t) is plotted as function 
of the time for different values of the monomer concen- 
tration. The data have been averaged over 40 different 
initial configurations. The different curves correspond to 
different values of cj), ranging from 0.5 to 0.85. At low 
concentrations the system completely relaxes within the 
simulation time; the relaxation process becomes slower 
as the concentration is increased and above the percola- 
tion threshold, </> c , the system is kinetically arrested, in 
the sense that the relaxation functions do not go to zero 
within the time scale of the simulations. 

We analyze more quantitatively the long time decay 
of fq(t) in Fig|3| As the monomer concentration, </>, ap- 
proaches the percolation threshold, <p c , fq*(t) displays a 
long time decay well fitted by a stretched exponential 
law - e-^l-rf with a ~ 0.30 ± 0.05. Intuitively, this 
behavior can be related to the cluster size distribution 
close to the gelation threshold, which produces relaxation 
processes taking place over different length scales. At 



the percolation threshold the onset of a power law decay 
(with an exponent c) is observed as shown by the double 
logarithmic plot of Fig|3 6]. As the monomer concen- 
tration is increased above the percolation threshold, i.e. 
in the gel phase, the long time power law decay of the 
relaxation functions can be fitted with a decreasing expo- 
nent c, varying from c ~ 1.0 at 4> c to c ~ 0.2 well above 
(f> c , where a nearly logarithmic decay appears. This sug- 
gests that the growth of the relaxation time is driven by 
the formation of the critical cluster, that actually deter- 
mines the kinetic arrest. On the whole, this behavior 
well agrees with the one observed in gelling systems in- 
vestigated in the experiments of refs. Q ■ It is interesting 
to notice that this kind of decay with a stretched expo- 
nential and a power law reminds the relaxation behavior 
found in spin-glasses |32j | . Many analogies in the dynam- 
ics of gels and spin-glasses have been recently pointed out 
[23, but the underlying physics is rather unclear. 

In the model with finite lifetime bonds, the equili- 
bration time is an increasing function of 77,. The sys- 
tem is considered at equilibrium when both the average 
number of bonds has reached its asymptotic limit and 
the autocorrelation functions have gone to zero [34|- In 
this case very different behaviors are observed. In FigQ] 
fq*(t) is plotted as function of time for a fixed value of 
n = 10, 100, 1000 (respectively FigEh,H> andgt) for in- 
creasing values of the monomer concentration (4> varies 
on the same range as the permanent bond case). At low 
concentrations, the behavior of the autocorrelation func- 
tion fq(t) is apparently very similar to the one observed in 
the case of permanent bonds: the system completely re- 
laxes within the simulation time scale and the relaxation 
time increases with the concentration cj). At high con- 
centrations, a two-step decay appears, closely resembling 
the one observed in super-cooled liquids. This qualita- 
tive behavior is observed for many different values of the 
bond lifetime, 77,. As shown in Fig. the two step decay 
is more pronounced for higher values of 17, . 

As we can see in Fig0 where the long time decay of 
fq(t) for 77, = 100 is shown, the long time decays are 
well fitted by stretched exponentials. The exponent 
(0 ~ 0.7 for the case of Fig|5j does not seem to vary sig- 
nificantly as the concentration varies, and this has been 
observed for all the values of 17, studied. It instead de- 
creases as 77, increases: for very small 77, one recovers a 
long time exponential decay whose behavior becomes less 
and less exponential as the bond lifetime increases. This 
suggests that the stretched exponential decay is due to 
the presence of long living structures. 

For high monomer concentrations we fit the fq(t) 
curves using the mode-coupling /3-correlator [T^. corre- 
sponding to a short time power law ~ / + C^J an d 

a long time von Schweidler law ~ / — (jn) ' ^ n Fig ED 
we show the agreement between the fit (the full lines) 
and the data for 77, = 1000 in the range of concentrations 
4> = 0.8 — 0.9. The exponents obtained are a ~ 0.33±0.01 
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and b ~ 0.65±0.01. At long times the different curves 
obtained for different <f> collapse onto a unique master 
curve by opportunely rescaling the time via a factor t(0) 
(Fig0 . The master curve is well fitted by a stretched ex- 
ponential decay with j3 ~ 0.50±0.06. The characteristic 
time t(4>) diverges at a value <\> g ~ 0.963±0.005 with the 
exponent 7 ~ 2.33±0.06 (Fig[SJ). This value well a gree s 
with the mode-coupling prediction 7 = l/2a+ 1/26 [l2j . 

The same behavior and the same level of agreement 
between the data and the mode coupling predictions 
have been obtained for different values of (77, = 
100, 200, 400, 1000, 3000) . Neither the exponents a and 
b obtained by the /3-correlator nor the exponent (3 of the 
stretched exponential vary significantly as function of <f> 
and of Tb. The value of <j> g , where the characteristic time 
r((f)) apparently diverges, seems instead to vary with Tfc, 
but it is always close to <f> g = 1. 

This glassy relaxation pattern suggests that also in 
this case the relaxation takes place by means of a caging 
mechanism: particles are trapped in a cage formed by 
the surrounding ones, the first relaxation step is due to 
movements within this cage, whereas the final relaxation 
is possible due to cage opening and rearrangement. We 
notice that, contrary to the usual behavior observed in 
super-cooled liquids and predicted by the Mode-Coupling 
Theory, the value of the plateau of the relaxation func- 
tions, which is typically related to the size of the cage, 
strongly depends on the monomer concentration, 4>. This 
effect will be explained later in terms of effective clusters. 



IV. THE RELAXATION TIME AND THE ROLE 
OF EFFECTIVE CLUSTERS 

We study now the relaxation times that can be ob- 
tained from the fq(t), as the time r such that fq{r) ~ 0.1. 
In Fig|U| the relaxation time r is plotted as function of 
the monomer concentration, <f>, for the permanent bonds 
and for the finite lifetime bonds case at different values 
of Tb- In the figure one finds the data for the permanent 
bond case on the left, and then from left to right the data 
for decreasing values of bond lifetime, Th- 
in the case of permanent bonds (chemical gelation), 
t{4>) displays a power law divergence at the percolation 
threshold <f> c . For finite bond lifetime the relaxation time 
instead increases following the permanent bond case, up 
to some value <jf and then deviates from it. The longer 
the bond lifetime the higher 0* is. For higher <f> the fur- 
ther increase of the relaxation time corresponds to the 
onset of the glassy regime in the relaxation behavior dis- 
cussed in the previous section. This truncated critical 
behavior followed by a glassy-like transition has been ac- 
tually detected in some colloidal systems in the viscosity 
behavior P, El- 

In both cases of permanent bonds and finite lifetime 
bonds, clusters of different sizes are present in the sys- 
tem. In the permanent bond case, a cluster of radius 
R diffuses in the medium formed by the other percola- 



tion clusters with a characteristic relaxation time t(R). 
At the percolation threshold the connectedness length 
critically grows in the system and so does the overall re- 
laxation time. In the case of a finite bond lifetime tj,, it 
will exist a cluster size R* so that Tb < t(R*). That is, 
clusters of size R> R* will break and lose their identity 
on a time scale shorter than t(R) and their full size will 
not contribute to the enhancement of the relaxation time 
in the system. We can say that the finite bond lifetime 
actually introduces an effective cluster size distribution 
with a cut-off and keeps the macroscopic viscosity finite 
in the system |37j . 

At high concentrations the system approaches a glassy 
regime and the relaxation time increases. In order to 
further investigate the high concentration regime, in 
FigEH we directly compare f g (t) at fixed <j> = 0.85 for 
71 = 10,100,400,1000, and the permanent bond case. 
We observe that at a fixed value of the monomer concen- 
tration, <f>, the relaxation functions calculated for finite 
lifetime bonds coincide with the permanent bond case up 
to times of the order of Tb- This suggests that on time 
scales smaller than n the relaxation process must be on 
the whole the same as in the case of permanent bonds, 
where permanent clusters are present in the system, and 
gives an interpretation in terms of effective clusters for 
the two step glassy behavior of the relaxation functions: 
The first step is due to the relaxation of a cluster within 
the cage formed by the other clusters, whereas the sec- 
ond step is due to the breaking of clusters. This second 
relaxation is the analog of the cage opening in an ordi- 
nary supercooled liquid. In conclusion, on a time scale 
of the order of Tb, the effective clusters play the role of 
single molecules in an ordinary supercooled liquid, or in 
a colloidal hard sphere system. 

Using this picture of effective clusters, we are able to 
explain the increase of the plateau in f q (t). In fact, since 
different values of the monomer concentration correspond 
to different effective cluster size distributions, for each 
value of 4> one has a different glassy liquid of effective 
clusters. This will change the first relaxation and should 
correspond to a change in the plateau of the relaxation 
functions (FigQ] and [fJJ . In particular we find that for 
higher <f> one has a higher plateau, that is the first de- 
cay (the motion of clusters within the cages) produces a 
smaller relaxation in the system. 



V. DIFFUSION PROPERTIES 

In order to obtain further information on the dynam- 
ics we calculate the mean square displacement of all the 

particles, (f*(t)) = j, E^i<(^(* + O " W)?)- ^ the 
model with finite lifetime bonds, clusters continuously 
evolve in time and therefore the diffusion coefficient of a 
single cluster cannot be defined. 

In the model with permanent bonds the mean square 
displacement of the particles (r 2 ^)) presents a long-time 
diffusive behavior, and the diffusion coefficient decreases 
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but remains finite also above <f> c . This is due to the fact 
that the infinite cluster can be viewed as a net with a 
large mesh size, through which monomers and small clus- 
ters can diffuse. 

In previous papers the diffusion coefficient of clus- 
ters with a fixed size was studied. We found that clusters, 
whose size is comparable with the connectedness length, 
present a diffusion coefficient going to zero at <pc (with 
the same exponent as the relaxation time), whereas single 
monomers present a finite diffusion coefficient also in the 
gel phase. As we have already noticed, this is due the fact 
that small clusters are able to escape through the perco- 
lating cluster having a structure with holes over many dif- 
ferent length scales close to the percolation threshold. It 
is therefore clear that the behavior here observed for the 
mean square displacement is mainly due to the diffusion 
of single monomers and small clusters. In Figllll {r(fij) 
is plotted as function of time in a double logarithmic plot 
for (j> approaching </> c in the case of permanent bonds: 
particles can still diffuse, and the diffusion coefficient ap- 
parently decreases with the monomer concentration. At 
high concentrations the sub-diffusive regime stays longer, 
and the long time diffusive behavior is hardly recovered. 

In Figll2l we plot the data obtained with r\, — 1000 
and <j) = 0.8, 0.82, 0.85, 0.9. According to the results just 
discussed, for low concentrations (r 2 ^)) shows a simple 
diffusive behavior. The diffusion does not change signifi- 
cantly close to 4> c and for high concentrations the behav- 
ior observed, characterized by a plateau, is similarly to 
glass forming systems. This onset of a glassy regime has 
been obtained for different values of Tb, and again it indi- 
cates a caging mechanism in the dynamics. The asymp- 
totic diffusion coefficient goes to zero as cj> approaches <p g 
(inset of Fig ll2ll . as a power law, with an exponent close 
to 7 (section ITV|) in agreement with the Mode Coupling 
Theory predictions. 

As already done for the relaxation functions we directly 
compare the mean square displacement obtained in the 
cases of permanent bonds and finite bond lifetime. Fig ll3l 
shows for a fixed value of the concentration, <j> = 0.85, 
that the two quantities coincide up to time scales of the 
order of t&. For longer times in the system with non 
permanent bonds the final diffusive regime is recovered. 
These results are coherent with the behavior of the re- 
laxation functions discussed in the previous section. The 
first regime is apparently related to the diffusion of effec- 
tive clusters. Here again the value of the plateau in the 
diffusion pattern, which is related to the size of the cage, 
varies with the concentration (f>. For higher values of the 
monomer concentration, the size of the cage apparently 
decreases. This corresponds to larger effective clusters, 
which have less free space compared to their size. By 
means of the qualitative argument used in section IIVI 



one expects that for a longer Tb the condition Tb < t(R*) 
will be fulfilled by a larger size R* , and on average larger 
clusters will persist. For the same value of the concen- 
tration, the size of the cage should be the same, whereas 
the particles of this glassy system (i.e. the effective clus- 
ters) are longer trapped in the cage as the bond lifetime 
increases fFig ll3|) . 



VI. DISCUSSION AND CONCLUSION 

We have studied a model for gelling systems both in 
the case of permanent bonds and finite bond lifetime. 
The study of the dynamics shows that when bonds are 
permanent (chemical gelation) the divergence of the re- 
laxation time is due to the formation of a macroscopic 
critical cluster and the decay of the relaxation functions 
is related to the relaxation of such cluster. In the case 
of finite Tb there is an effective cluster size distribution, 
with a size cutoff. Note that the clusters cannot be easily 
defined as in the case of chemical gelation: The effective 
clusters do not coincide with pairwise bonded particles. 
A cluster can be identified in a statistical sense as a group 
of monomers which keeps its identity (i.e. the bonds 
are unbroken) when diffusing a distance of the order of 
its diameter. The formation of effective clusters leads 
to an apparent divergence of the relaxation time which 
is characterized by exponents corresponding to the case 
of random permanent bonds (random percolation). As 
the monomer density increases the presence of effective 
clusters further slows down the dynamics, until a glass 
transition is reached. 

In the case Tb — > oo, starting with a random config- 
uration of unbonded monomers one obtains an out of 
equilibrium state as in cluster-cluster aggregation mod- 
els, which can be representative of colloidal gelation. Ide- 
ally this out of equilibrium system is connected to the two 
lines described above, the pseudo percolation line and the 
glassy line. The pseudo percolation line can be detected 
if the effective cluster size is large enough and it is quite 
distinct from the glassy line. However both lines inter- 
fere at low densities and low temperatures with the phase 
coexistence curve. 
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FIG. 1: Two typical configurations obtained for <f> = 0.6 
with Tb — 100 and /(, = 0.02 (a), where there is no evidence 
of phase separation, and with r b — 2, f b — 1 and monomers 
of valence 6 (6), where the phase separation occurs. 



FIG. 2: The relaxation functions for permanent bonds fq(t) 
as function of the time for q ~ 1.36 and, from left to right, 
(j> = 0.5, 0.6, 0.68, 0.718, 0.75, 0.8, 0.85. 
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FIG. 3: Double logarithmic plot of the autocorrelation 
functions fq(t) as function of the time for q ~ 1.36 and 
<t> = 0.6, 0.718, 0.8, 0.87. For cj> < <f> c the long time decay is 
well fitted by a function (full line) ~ e" (t/T) ' 3 with /3 ~ 0.3. 
At the percolation threshold and in the gel phase in the long 
time decay the data are well fitted by a function ~ ( 1 + -p- ) ~ c . 
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FIG. 4: fq(t) as function of the time for q ~ 1.36 
calculated on a cubic lattice of size L = 16: for <j> — 
0.6,0.7,0.8,0.85,0.87,0.9 (from left to right) and n = 
10MC step /particle (a); Tb = 100MC 'step /particle (b); n = 
1000MC step /particle (c); the dotted lines are a guide to the 
eye. 




FIG. 5: The long time decay of f$(t) in a log-log plot for 
q ~ 1.36. It has been calculated on a cubic lattice of size 
L = 16 for r b — 100 MC step /particle (from left to right 
4> = 0.65,0.68,0.718). The data are fitted using a stretched 
exponential function with f3 ~ 0.7 (full lines). 
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FIG. 6: fq{t) as function of the time for q ~ 
1.36 calculated on a cubic lattice of size L = 16 for 
r b — WOOMCstep/particle. From left to right <j> = 
0.8, 0.85, 0.87, 0.9. The full lines correspond to the fit with 
the /3-correlator. 




FIG. 7: fq(t) obtained for q ~ 1.36, r b = 1000, and 4> = 
0.85, 0.87, 0.9, 0.91, 0.92: by opportunely rescaling them by a 
quantity r(0) they collapse into a unique master curve, well 
fitted by a stretched exponential function with f3 ~ 0.5 ±0.06. 
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FIG. 8: Log-log plot of the characteristic time t{4>) obtained 
by the rescaling of the relaxation functions. The points have 
been fitted (full line) with the function 0.006(0,, - 0)~ 2 ' 33 , 
where <p g ~ 0.96 ± 0.01 . 
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FIG. 9: The average relaxation time as function of the 
density; from left to right: the data for the permanent 
bond case diverge at the percolation threshold with a power 
law (the full line); the other data refer to finite ti = 
3000, 1000, 400, 1QQMC step /particle decreasing from left to 
right (the dotted lines are a guide to the eye). The apparent 
divergence of the relaxation time, r, is observed at the perco- 
lation threshold of the permanent bond case, <j) c = 0.718, for 
all the values of r b . 
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FIG. 10: fq(t) obtained for 4> = 0.85 and q ~ 1.36: the 
different curves refer to Tt = 10, 100, 400, 1000, compared to 
the permanent bond case (from left to right). 
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FIG. 11: The mean-square displacement (f^(t)) of the par- 
ticles as function of the time in a double logarithmic plot for 
permanent bonds: from top to bottom <j) = 0.4, 0.5, 0.6, 0.7, 
approaching <j> c . 
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FIG. 12: The mean-square displacement of the par- 

ticles as function of the time in a double logarithmic plot 
for t;, = 1000MC step /particle: from top to bottom (j> = 
0.8, 0.82, 0.85, 0.87, 0.9, approaching (f> g (T b ). In the inset, the 
diffusion coefficient: the full line is the fit with the function 
~ (0.963 - (j>)-' 2 - 3 . 
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FIG. 13: The mean-square displacement (r°(t)) of the par- 
ticles as function of the time in a double logarithmic plot, 
obtained at <f> = 0.85. The different curves, from top to bot- 
tom, refer to r& = 10, 100, 1000 and the case of permanent 
bonds. 



